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Abstract 

A popular observable in finite-temperature lattice QCD is the so-called singlet quark-anti- 
quark free energy, conventionally denned in Coulomb gauge. In an effort to interpret the 
existing numerical data on this observable, we compute it at order 0{a 2 s ) in continuum, and 
analyze the result at various distance scales. At short distances (r <C 1/nT) the behaviour 
matches that of the gauge-independent zero-temperature potential; on the other hand at large 
distances (r S> 1/ vrT) the singlet free energy appears to have a gauge- fixing related power-law 
tail. At infinite distance the result again becomes physical in the sense that it goes over to 
a gauge-independent disconnected contribution, the square of the expectation value of the 
trace of the Polyakov loop; we recompute this quantity at 0{a 2 s ), finding for pure SU(iV c ) a 
different non-logarithmic term than in previous literature, and adding for full QCD the quark 
contribution. We also discuss the value of the singlet free energy in a general covariant gauge, 
as well as the behaviour of the cyclic Wilson loop that is obtained if the singlet free energy 
is made gauge-independent by inserting straight spacelike Wilson lines into the observable. 
Comparisons with lattice data are carried out where possible. 
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1. Introduction 



One of the classic probes for forming a thermalized partonic medium in a heavy ion collision is 
the change that this should cause in the properties of heavy quarkonium pQ . It is traditional 
to address quarkonium through potential models, but in the finite-temperature context this 
has been hampered by the multitude of independent definitions of a potential that could in 
principle be introduced. A possible way to reduce the degeneracy is to use the weak-coupling 
expansion as a test bench; at least within that approach, it should be possible to derive from 
QCD the potential that plays a physical role. Indeed, it has been discovered that the relevant 
potential is none of the plethora of conventional ones but a new object, which even has an 
imaginary part, representing a decoherence of the quark-antiquark state that is induced by 
the thermal medium [2] _ [5]. Subsequent to the perturbative introduction, it appears that it 
may be possible to promote the corresponding definition to the non-perturbative level [6]. 

In this paper, we do not discuss the question of a proper definition of a static potential at 
finite temperatures, but rather set a more modest goal. Extensive lattice measurements have 
already been carried out on many potentials (see, e.g., refs. [7]-|15j and references therein), 
and first tests [6] suggest that the real part of the proper potential may overlap with one of the 
existing objects, namely the "singlet" quark-antiquark free energy. This concept is, however, 
based on gauge fixing, more precisely on the Coulomb gauge (cf. eq. (|3.ip below). We are not 
aware of a satisfactory a priori theoretical justification for this choice; the original motivation 
appeared rather to be practical, in that the Coulomb gauge potential empirically reproduces 
the expected zero-temperature behaviour at short distances (see, e.g., refs. [HUE]), and also 
displays manageable statistical fluctuations as well as a good scaling with the volume and 
the lattice spacing [8]. 

Given these empirical observations, we feel that there may be room also for some theoretical 
work of the Coulomb gauge potential. The purpose of this paper is straightforward: we 
compute the Coulomb gauge potential, and several variants thereof, at next-to-leading order 
in the weak-coupling expansion (i.e. 0(g 4 ), where g 2 = Aira s is the QCD coupling constant) 
in dimensional regularization. The results turn out, indeed, to yield a qualitative surprise 
that could have a bearing on the interpretation of the Coulomb gauge lattice data. 

The paper is organized as follows. In sec. [2] we compute the expectation value of a single 
Polyakov loop in dimensional regularization, and compare the result with literature. In 
sec. [3] we analyze the Coulomb gauge singlet quark-antiquark free energy at various distance 
scales, identifying both fortunate physical aspects as well as what looks like an unfortunate, if 
numerically small, gauge artifact in the result. In sec. U] we briefly discuss the corresponding 
object in a general covariant gauge, while sec. [5] looks into the so-called cyclic Wilson loop; 
this gauge-invariant completion removes the gauge artifact from the Coulomb gauge result 
but with the price of introducing new problems. Sec. [6] summarizes our main findings. 
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2. Polyakov loop expectation value 



2.1. Basic setup 

Employing conventions where the covariant derivative reads D„ = — ig B A^, with g B the 
bare gauge coupling and a traceless and hermitean gauge field, we define a "Polyakov 
loop" at spatial position r through 

P r = 1 + igs f drA (r, r) + {ig B f f 'dr f dr' A {r, r)A (r', r) + . . . . (2.1) 
Jo Jo Jo 

It transforms in a gauge transformation U(r,r) as P r — > U(j3, r) P r C/ _1 (0, r). Likewise, for 
future reference, a straight spacelike Wilson line from origin to r, at a fixed time coordinate 
r, is denoted by 

rl rl rX 

W T = t + ig B d\r-A(T,\r) + (ig B ) 2 dX / dA' r • A(r, Ar) r • A(r, A'r) + . . . , (2.2) 
Jo Jo Jo 

and transforms as W T — > U{t,y)W t U~ 1 {t,0). Note that because of the usual periodic 
boundary conditions in the time direction, Wp = Wq. 

In order to simplify the notation somewhat, we in general leave out the subscript from 
g B in the following formulae. It is to be understood that at the initial stages it is the bare 
coupling that is referred to, while at the end the bare coupling is re-expressed in terms of the 
renormalized one, 

2 2 g y- 2e 2Nj - UN C 6A 
9b=9 + -j^yr j e + 0(g), (2.3) 

where g 2 denotes the renormalized gauge coupling in the MS scheme, with a scale parameter 
jx 2 = 47re~ 7E /U 2 . The factor pT 2<1 is normally not displayed explicitly. 

Within perturbation theory, assuming that the Z{N C ) center symmetry is broken in the 
trivial direction, either spontaneously or explicitly due to matter fields in the fundamental 
representation, we can then define the expectation value of the trace of the Polyakov loop 
through [16J 

V^p = -^(Tr[P r ]> , (2.4) 



N c 

where the brackets refer to averaging in a thermal ensemble at a temperature T. On the 
non-perturbative level, ip P plays a role as the disconnected part of some 2-point function, 
e.gEl 

l^l 2 = , ! im ^2 ( Tr Tr [P f ]) = Urn -^{Tr [P r P f ]) coulomb • (2.5) 
r|^-oo iVjf |r|— s>oo iv c 



1 The second equality can be understood by writing P T = J- Tr [P T ] + P r , where P T is traceless and changes 
in gauge transformations like an adjoint scalar; P r should not have correlations at infinite distance, and our 
results confirm this expectation at 0{a 2 s ). 
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Figure 1: The graphs contributing to the expectation value of a single Polyakov loop up to 0(<? 4 ), 
with the filled blob denoting the 1-loop gauge field self-energy. 



Any practical computation in finite-temperature QCD is hampered by infrared divergences. 
There are two kinds: milder ones, associated with colour-electric modes and the scale gT, 
which can be cured by a certain resummation of the perturbative series; and more serious ones, 
associated with colour-magnetic modes and the scale g 2 T/n, which can only be handled by a 
non-perturbative study of three-dimensional pure Yang-Mills theory [17]. For the observables 
that we are concerned with, the first type plays a more prominent role; it can be handled by 
writing the observable as 



pJqcd 



p;eqcd 



+ 



unrcsummcd 



EQCD 



(2.6) 



The difference in the first square brackets is infrared finite (provided that the correct low 
energy effective theory [18], called EQCD [15] ; is used), and can be computed in naive per- 
turbation theory. The second term is infrared sensitive, but the resummations required for 
computing it can be implemented in the simple three-dimensional EQCD framework. 

The graphs contributing to the expectation value of a single Polyakov loop up to 0(g 4 ) are 
shown in fig. [TJ The computation of the graphs is in principle straightforward, however some 
care is needed when treating Matsubara zero modes. We separate their effects explicitly, 
writing 



drV^' = * fc T+(l-fJ; 



I _/ Jq n r 



dr' t' e 



iq n 



§ qn y + (1 - Sq n ) 



T 

iq n 



e iqnT + — I < 



iqnT 



(2.7) 
(2.8) 



etc, where q n is a bosonic Matsubara mode, q n = 2nTn, and 5 Qn = 5 n> o is a Kronecker delta 
function. Adding up the graphs in any gauge where the gauge field propagator, G^ v , has the 
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property Goi(0, k) = 0, the result (for the moment unresummed) can be written aj§ 
ip P = exp 



1 



l - g 2 C F p ^Goo(0, k) + 1 g 2 C F p Jc^(0, k) n 00 (0, k) 
/3 2 G 00 (0,k)Goo(0,q) 



■g*C F N ( 



k.q 



24 



■^[-Goo(0,k) + iGoo(«n,k) 



-2^G 00 (0,k) 



Goo(<?n,q) 



Goo(gn,q) 
1 + G i(q n ,q) 



G 0i ((?„,q + k) I , (2.9) 



where C F = (iV^ — l)/2iV c ; (3 = 1/T; H^ v is the gauge field self-energy; and J2q> = J2q n ^o- 
Inserting the propagator of the general covariant gauge (eq. (|4.2p below) as well as the 
corresponding self-energy (eq. (14. 3p below), it can be checked that the value of the expression 
in eq. (12. 9\\ is gauge-parameter independent. The same result is also obtained in the Coulomb 
gauge, by making use of eqs. (j3.3l) . (13.40 . In explicit form, the gauge-independent part can 
be written as 



■0F 



g 2 C F (3 f 1 1 fg 2 C F p f 1 



k 2 + 2 



k 2 



g A C F N t 



g 4 C F N c 



E 

{in} 

£ 



+ 



+ 



k, Q L W k*Q 2 (Q + K) 2 k 2 Q 2 (Q + K) 2 

~2_-D_ 2(D - 2)q 2 n 2 
k, q L W k*Q 2 {Q + K) 2 ^ k 2 Q 2 (Q + Kf 

1 1 



g 4 C F N c 



k,q 



2^Q2 (? 2 +A .2) fc 2(? 2Q2 

/? 2 



2 — D 



+ 



k i q 2 k 2 q 2( q + k y 2Ak 2 q 2 



(2.10) 



where = |k|, q = |q|, fc + q = |k + q|; K = (0,k), Q = (q n ,<i)', and we have separated the 
Matsubara zero mode contribution in the term of 0(g 4 ). 



2.2. Soft-mode contribution 

It is immediately visible from eq. (|2.10p that the k-integral in the term of 0(g ) is infrared 
divergent, ~ J^l/k 4 . The coefficient of the divergence is, however, nothing but the Debye 

2 Here and in the following we often use exp(x) as a shorthand for 1 + x + x 2 /2; no proofs concerning the 
exponentiation of higher order corrections are provided. 
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mass parameter: 



g 2 c F /3 r 1 



... + 



{q n } 

ml + 0(k 2 ) 



2 A k 1 

where we made use of the sum- integrals in eqs. (EOl) . (HU) , and denoted 



6 3 / 



(2.11) 



(2.12) 



Therefore the divergence can be removed by the usual resummation of colour-electric modes, 
g 2 C F f3 f 1 , g 2 C F p f m 2 E , 5 2 C F /3 /" 1 g 2 C F m E (3 



1 ff 2 CF/3 /-ml 



+ . . . 



k A; 2 + m 



8tt 



(2.13) 



Unfortunately, an inspection of the last row in eq. ()2.10|) shows that there are also loga- 
rithmic infrared divergences (cf. ref. |20|). and to handle them correctly we need to proceed 
carefully. A systematic way is through eq. ()2.6|) . Within EQCD, the effective Lagrangian has 
the form 

(2.14) 



\ Tr [F 2 ] + Tr [D U A ] 2 + m 2 Tr [A 2 ] + ... , 



where F {j = (i/g E ){Lh, Dj], L>i = d { - ig E Ai, k\ = A°;T a , A = A%T a , and T a are hermitean 
generators of SU(3)c| The Polyakov loop operator is represented as 



igA p Z 1 + - (igAoP) 2 Z 2 + . . . + (g 2 ^? 2 ) 2 ^ + ... 



(2.15) 



Here all possible local operators made of Matsubara zero-modes, invariant under parity and 
spatial rotations and transforming under the adjoint representation of the gauge group, can 
in principle appear. It will be convenient for our purposes to use a "mixed" convention in 
eq. ()2.15p where g denotes the (renormalized) four-dimensional coupling while the fields are 
those of three-dimensional EQCD. The matching coefficients Z{ are of the form Zi = l+0{g 2 ). 
The possible appearance of a matching coefficient like was pointed out in ref. |22j . but 
it does not contribute at the order of our computation. The first term in eq. (|2,15p has 
been put inside brackets, because in the language of eq. (12.6P it represents the value of the 
unresummed difference inside the first brackets of eq. (I2.6P . 

Now, within EQCD, all dynamical effects involve the scale m E ~ gT and thus bring in 
additional powers of the coupling g. In fact, the leading term, originating from the 2nd order 
operator in eq. (I2.15p . precisely reproduces the result of eq. (I2.13P which is of order 0(g 3 ). 



3 A 2-loop derivation of m%, gj| in terms of the parameters of four-dimensional QCD can be found in ref. [21] , 
but here we only need their values at leading non-trivial order, m| = g 2 T 2 (Nf/6 + 2V c /3), g\ = g 2 . 
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Because the term is of 0(g 3 ), the coefficient Z2 can be set to unity as we work at the order 
0(g 4 ). There is an 0(g 4 ) contribution from the next-to-leading order evaluation of the 2nd 
order operator, however, and the sum of the 0(g 3 ) and 0(g 4 ) contributions can be written 

as 



(V>p) 



EQCD 



g 2 C F m E p , g 2 C F p 



+ 



8vr 2 h(k 2 +ml) 2 

where the self-energy within EQCD has a well-known form (see, e.g., refs. [231 12H E5]) : 

1 



g 2 E N c T 



+ 



2(m£ - k 



q 2 + m 2 g 2 [(/c + g) 2 + m|] 



+ 



+ - 



D 



+ 



£(k 2 + 



mi 



k + m 



(k + g) 2 + m : 



2.16) 



(2.17) 



We have kept a general gauge parameter £ here, defined with the convention that £ = — 1 
corresponds to Coulomb gauge and £ = to Feynman gauge (i.e. G{j = 5ij/q 2 + S, q%qj /q )', 
and, for later reference, we have not yet killed any integrals through special properties of 
dimensional regularization. 

We note, first of all, that the ^-dependent part of eq. (|2.17p gives no contribution in 
eq. (|2.16|) . Second, the other terms of eq. (|2.17p reproduce, for m E — > 0, the zero- mode 
contribution on the last line of eq. (|2.10p . In the "dressed" form of eqs. (|2.16p . (|2.17p . 
however, the logarithmic infrared divergence has been lifted. The integrals in eq. (|2.16p ar e 
all elementary (cf. eqs. (|A.26|) ~ (jA.28|) ) . and in total we get (replacing g 2 — > g 2 + 0(g A )) 



P^EQCD 



g 2 C F mv(3 g A C F N c ( 1 



rcsummcd 



8tt 



(47T) 



4e 2m E 



+ \ ) + 0{g^ 



(2.18) 



2.3. Hard-mode contribution 



It remains to compute the unresummed difference inside the first brackets in eq. (|2.6p . given 
by eq. (|2.1U|) . Without resummation the zero-mode contribution in eq. (|2.1U|) contains no 
scale and vanishes; the same happens in the cases where the k-integration contains no q 2 , and 
factorizes from the q-integration. This leaves us with 5 non-zero sum- integrals, with values 
given in eqs. ()A.8p - (jA.12p ; summing together, we arrive at 



Za 



(V>p) 



1 + 



QCD 



9 A C F 



(V>p) 



EQCD 



unrcsum mod 



„,_-) + , (J_ + ln A4 



(2.19) 
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2.4. Summary and comparison with literature 

Adding up eqs. (12. 18ft . (|2.19j) . the 1/e's duly cancel, and we get our final result for the 
Polyakov loop expectation value: 



g 2 C F m E , g 4 C F 



8vrT (4vr) 



, ln2\ / m B 1 



+ 0(9 5 ), (2.20) 



where m E is from eq. f)2 . 12|) □ In accordance with general expectations |26[ [57] , the result is 
finite and renormalization group invariant up to the order computed^ 

A classic determination of ip P was presented, for Af = 0, in ref. [20 . Rc-cxprcssing that 
result in terms of m E it can be written as 



, m 1 , g 2 C F m E ( g 4 C F N c m E i 3 



^ = 1 + » ~ r + g ~ \ o In — + - . (2.21) 
^ P 8vrT (4vr) 2 \ 2T 4 / v ; 

This agrees with our result including the logarithmic term, but differs on the constant term 
accompanying the logarithm by a factor — In 2 + 1/2. As far as we can see, the difference can 
be traced back to the way that the resummation was carried out. We do believe ours to be 
a systematic resummation. 

In order to finally compare with lattice results, we need to insert some numerical values 
for the parameters appearing in eq. (|2.20p . Following ref. [29] , we can estimate them through 
some "fastest apparent convergence" criteria; for the gauge coupling we have applied this 
criterion to the combination g 2 Z 2 , playing a role in the next section (cf. eqs. (|3.19p . (|3.28p ). 
while for the mass parameter we take over the criterion from EQCD [29J: 



9 ^ 2 ±1 m 2^ ^(2N c + Ni )T' 

(llA c -2A f )[lng- 7B + c g ] ' E (llA c -2A^)[lng- 7E + c m ] ' 

where 



2N t (4 In 2 - 1) - 11JV C 4A f In 2 5A C 2 + A f 2 + 9A f /(2A c 



(2.23) 



y 2(llA c -2A f ) ' HA c -2Af (11A C - 2A f )(2A c + N { ) ' 

The results are shown in fig. [2J where they are also compared with the four-dimensional 
lattice data from ref. [12J. We have also tested more elaborate choices for the parameters, 
leading indeed to a somewhat better accord with lattice data (cf. caption of fig. [2]), however 
within the accuracy of our actual computation it is not possible to justify these theoretically. 
Nevertheless experience from other quantities, such as the spatial string tension \\.\\ [21~1 [3TJ 
[32] . leads us to suspect that pursuing the computation systematically to a higher order would 
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If the quarks are given a common chemical potential, /x, then the Debye mass parameter gets changed as 
! [iVc^ + N f + ^t)], and the numerical factor in eq. (f2T20]i is modified as - In 2/2 -> Re C'(0, \ + 



*2^f)' wnere C( x i v) = 9 x (^(x, y), and ((x, y) denotes the generalized zeta function, (^(x, y) = J^feLo V(^ + 2/) a 
5 We have been informed by the authors of ref. [28] that they have recently obtained the same result. 
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Figure 2: The dimensionally regularized Polyakov loop expectation value, eq. (|2.20p . as a function 
of T/T c , where T c is the critical temperature of the deconfining phase transition (a conversion from 
perturbative units has been carried out by assuming T c /A^jg ~ 1.25; a variation within the range 1.10 
- 1.35 yields an error much smaller than the band width). The upper edge of the band (solid red line) 
corresponds to evaluating the coupling and the Debye mass parameter according to the simple 1-loop 
criteria in eqs. (|2.22[) . (|2.23j) : for the lower edge (dashed blue line) we have replaced g 2 through the 
2-loop value of g\ given in ref. [21], and m^/g 2 , through the expression in eq. (14) of ref. [30]. The 
lattice data, labelled by N T , is from ref. [H] (the spatial lattice size was kept fixed at 32 3 ). 

eventually allow to improve on the agreement. (Phenomenological recipes for matching the 
lattice data down to lower temperatures can be found, e.g., in refs. |33j.) 

It is amusing to note, in any case, that the behaviour of the Polyakov loop is qualitatively 
quite similar to that of mesonic screening masses, expressed in units of the temperature |34j : 
both are small close to the phase transition (because they are related to order parameters in 
various limits), but then increase rapidly, and should finally approach their non-zero asymp- 
totic values from above. 

3. Singlet free energy in Coulomb gauge 
3.1. Basic setup 

The original definition of the singlet quark-antiquark free energy was related to the eigenvalues 
of the untraced Polyakov loop [35]; in practice, however, lattice measurements consider the 



8 



VWWVWVA 



wvwww 



wwwwwwv 



WVWWWV 



WAAAMM 



1 ^Jjlfww 



Figure 3: The connected graphs contributing to the correlator of two Polyakov loops at 0(<? 4 ), 
with the filled blob denoting the 1-loop gauge field self-energy. In addition there are disconnected 
contributions, obtained by multiplying the set of graphs in fig. [TJ with its Hermitean conjugate. (The 
topologies are the same as those needed for pushing the computation of the heavy quark medium 
polarization one level up from the current standard [37].) 



object (see, e.g., [8]-[T5] and references therein) 



Mr) = ^r{Tr[P T Pi]) 



Coulomb 



(3.1) 



This differs from the physical "colour-averaged" free energy, in which traces of Polyakov loops 
are correlated (cf. eq. (|5.30|) ) [T6J, [231 ES] ; the advantage of eq. (|3. 1 j) is that the free energy 
extracted from it is believed to reduce at small distances to the gauge-fixing independent 
zero-temperature static potential [THUS]. 

The general strategy for our determination of ipc is the same as in sec. [2] The new graphs 
to be evaluated are shown in fig. [3] Like in the previous section, we first give a general 
(unresummed) result, which is valid in any gauge with the property Go«(0,k) = 0: 



ipc(r) = exp 



tfCpP /(e ikr -l)G o(0,k 



-g z C F p J[e %kr - 1) G z m (0, k) n 00 (0, k) 



qn^O 



G o{q n , q)— — — + G oi (q n , q) 



Qn 



G 0i (q n , c[ + k) 



/3 2 G 00 (0,q)Goo(0,k) 



+ E 



Goo(q n ,q) 



qn^O 



Goo(0,k) + -Goo(^,k) 



(3.2) 



This was obtained by a brute force evaluation, making use of eqs. (|2.7p . (|2.8p . etc. The 
result is gauge-dependent; in order to now specialize to Coulomb gauge, we insert the gluon 
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propagator 



■fl^l^O' r)' <2 = (<?«, q) 



Q 2 



(3.3) 



and the corresponding self-energy into eq. f|3.2|) . The self-energy reads 

L>-2 2[k 2 + (D - 2)g 2 ] 



n 00 (o,k) 



s 2 ^ 



A: 2 



+ 



Q 2 Q 2 {Q + K) 2 

k 4 k 4 q 2 



q 2 Q 2 q 2 (q + k) 2 Q 2 2q 2 (q + k) 2 Q 2 (Q + K) 2 



{Q} 



k 2 + 4g 2 
Q 2 ' Q 2 (Q + K) 2 



+ 



(3.4) 



where again K = (0,k). To be explicit, the result can be written as 

3 ik-p _ 1 1 



l + 9 z C F /3 



A; 2 



e ik-r _j\2 



A; 2 



+g^C F N c Y J l (e 4k ' r -l 



1 



4<Z 2 



A; 4 Q 2 k 2 Q 2 {Q + K) 2 k*Q 2 (Q + K) 2 



D 



+ 



+ 



2(D - 2)g 2 



A; 4 Q 2 k 2 Q 2 {Q + K) 2 k*Q 2 (Q + K) 2 
1 1 



k 2 q 2 Q 2 2q 2 Q 2 (q + k) 2 2q 2 Q 2 (Q + K) 2 



+g 4 C F N c 
+g 4 C F N c 



e ik-r _ 1 



A; 2 



6 



i 



gj(k+q)-r _ j 



k,q 



A; 2 



8 ^Z^y'g 2 



(3.5) 



Given that ^2 q , q~ 2 = /3 2 £(2)/27r 2 = (3 2 /l2, the last two terms (corresponding to a vertex 
correction and to two-gluon exchange) actually do not contribute; in Coulomb gauge the 
unresummed result originates from the self-energy correction alone. 

If we naively set |r| — > oo in eq. (|3.5p . then all terms containing e jkr or e i ( k+q )' r drop out. 
A constant term remains over, and as can be verified by comparing with eq. (|2.10p . it has 
the valuj§ 

lim Vc(r) = |V>p| 2 • (3-6) 

r— >oo 

In the following we drop out this constant contribution, and focus on the r-dependent terms. 



6 The expression in eq. (|2.10|) had the covariant gauge as a starting point whereas eq. (|3 . 5 1) comes from the 
Coulomb gauge; therefore some of the terms do not immediately look alike. On closer inspection, however, 
their difference integrates to zero in the r-independent term. 
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3.2. Short-distance limit 



The Matsubara sums in eq. (|3.5h can be transformed, in a standard way, to an integral 
representation, from which a zero-temperature part and a thermal part can be identified. 
The thermal part contains the scale T inside Bose-Einstein and Fermi-Dirac distribution 
functions and, for large momenta, is exponentially suppressed. Therefore we expect that at 
small distances, i.e. r -C Jp, the result can be obtained by simply replacing the Matsubara 
sums with the corresponding zero-temperature integrals, $ Q ,^ {g} — > Jq. (Of course, this 
expectation can be crosschecked later on from the more general result.) 

The advantage of the zero-temperature limit is that then all integrals can be carried out 
analytically. For the "covariant" structures in eq. (|3.5p we get, as usual, 



r q 2 q 

J q Q 2 (Q + k) 2 
i 

1 

q Q 2 (Q~+Kf 



l 



D 
0. 



1 



A: 2 



2Q 2 AQ 2 {Q + K) 2 



1 



(47T) 



7 + >4 + 2 



(3.7) 
(3.8) 
(3.9) 



where we once again inserted K = (0,k). As far as the "non-covariant" terms are concerned, 
it is convenient to group them as 



1 



Q q 2 Q 2 
1 

To 2 



i 

¥ 
i 

k 2 



l 



i 



S7T 



2 k 2 



- +ln— - 

e Ak 2 



2(q + k) 2 
1 

(q + k) 2 

2 



+ 



2(Q + K) 2 
1 



Q 



2q 2 Q 2 



(q + k) 2 {Q + K) 2 \ 



+ 



8ir 2 k 2 



In 4 . 



(3.10) 



The former combination is infrared finite but ultraviolet divergent, and the integral was 
carried out in dimensional regularization; the latter combination is infrared and ultraviolet 
finite for k ^ 0, and was integrated directly in four dimensions. 

Summing up and re-expanding the bare gauge coupling of the leading order term in terms 
of the renormalized one, the 1/e's cancel (more details on the graph-by-graph origin of diver- 
gences are provided in sec. H]for a case where they do not cancel), and we get 

r7rT<l ( g 2 CF/3\ 



exp 



47rr 



9 4 C F p fe 
(4tt) 2 



A: 2 



2N f 



1 P + 5 
ln ^ + 3 



+ 



llN r 



1 A 

ln F + 33 



(3.11) 



Exponentiating the 0(g 4 ) correction and extracting the singlet free energy, V±, as ipc( r ) = 
exp[—/3Vi(r)], we see that V\ precisely agrees with the classic result for the zero-temperature 
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static potential at 0(g 4 ) \38^\ : 



( \ 9 2 C F g 4 C F f e ikr [2iV f / fi 2 ,5\ UN C ( p, 2 31 \1 

Fi(r) = -1^7 + y k if Lt r ¥ + 3 J - — r f + 33 J J • (3 - i2) 

Empirical observations in the same direction have been made even on the non-perturbative 
level [141115]; thus the agreement might be true at the 2-loop level as well, where a comparison 
could be carried out with ref. [39], however we have not undertaken this task. 



3.3. Hard-mode contribution 

We now move on to consider the behaviour of eq. (|3.5p at larger distances, r ~ 1/vrT, where 
thermal effects do play a role. The general strategy is analogous to eq. (|2.6p . viz. 



(V'c(O)qcd - (V>c(0)eqcd + Oc(j))eqcd ■ (3.13) 



unrcsummcd 



rcsummcd 



This time we evaluate the unresummed QCD part first, and for ease of later subtraction 
separate the Matsubara zero mode contribution from eq. (|3.5p . i.e. replace Ylq n ~* > an< ^ 
write the zero-mode part separately. 

For the non-zero mode part, we proceed as follows. The goal is to carry out the sum- 
integral over Q of non-zero Matsubara frequency, and to express the resulting k-integrand, 
let us call it X(k), in the form 

l(k) = ^ + ^ + C(k). (3.14) 

Because of restriction to non-zero modes, the function C(k) must be analytic in k 2 . Like 
in eq. (|2.13p . the coefficient A is essentially m|, and this most infrared sensitive term is 
subtracted by the expanded version of the EQCD contribution. The coefficient B, in turn, 
"renormalizes" the leading order contribution; in fact it also gets subtracted by the EQCD 
contribution through eq. (|3.13p . more precisely by the 1st order term from eq. (|2.15p . with 
a properly chosen factor Z%. The remainder, determined by the function C(k), represents 
the "genuine" thermal contribution from the non-zero Matsubara modes, and does not get 
subtracted by EQCD effects. 



7 We recall that in the Fourier transform the logarithm gets effectively replaced as \n(fi 2 /k 2 ) — > 2[ln(/w)+7E]. 
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Making use of the sum-integrals in eqs. (fO|> . fO|) . (|A~7Tj) . (fQ3|) - ([Xl8|) . we get 



(^c(r)) 



QCD 



unresummed 
.2/ 



n,nst " 1 rA ' 3 / k V + 2 (/ Cf/? y k V 



+ 



k 4 (4vr) 2 A; 2 V 3 



^ (Lb + 1) _^£ (L[ _ i: 



+ 



+ 



9 A C F N C 
(4tt) 2 



9 4 C F N { 
(4vr) 2 



1 5 ln(l - e" 4 ^ r ) 1 + e 47ri r (4vrTr - 1 



24T 2 r 2 



QnTr 



+ ^(4- ( 47r ^ rn ) )/ ; !j7 ;/v,;) + 



n=l 



3 ( e 4^Tr _ 1)2 

Li 2 (e" 4 ^) 



1 



In 



3irTr 1 + e 



_ e -27rTr e 2irTr 
+ — =— 



-27rTY 



1 



(2ttTY) 2 

„4?rTr 



2vrTr 



( e 4^Tr _ !)2 



n=l 



[2ttTV x (2n - 1)1 2 , , 
-2 + ± ^ ^- ) J5 X 27rTr x (2n - T 



+g 4 CpN c x (zero mode contribution) . 



(3.15) 



Here the gauge coupling is already the renormalized one; Lb, Lf are defined as 



L b = 2 In 



7E 



4vrT 



, L f = 2 In 



fie 



7E 



vrT 



(3.16) 



and Ei and L12 are defined in eqs. (1X201) . (COD . Sometimes it may be convenient to 
replace the sums over Ei with integral representations, and indeed this can be achieved as 
specified in eqs. (|A.22|) - (|A.25|) . The expression then compactifies quite a bit; we rewrite the 
corresponding result in eq. (|3.28p below. 

As a first check, it can be shown that for irTr <C 1 the square-bracketed terms in eq. ()3. 15[) 
go over into 



g 4 C F N c 11 
(4vr) 3 Tr X T 



21n(4vrTr) - 1 + 



31 

33 



9 A C F N { 2 
(4vr) 3 Tr X 3 



21n(7rTr) + 1 + 



(3.17) 



Combining with the logarithmic terms in eq. (|3.15p . we exactly match the behaviour given 
by eq. (|3Tljh 
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3.4. Soft-mode contribution 



We then proceed to consider the EQCD (zero- mode) contribution to eq. (13,13p . which we 
now denote by ipc(r) = (tpc( r ))EQcu- It reads 



Mr) 



const. + g 2 C F pZ 2 



„ik-r 



k 2 + m 



1 

2+ 2 



+ g 2 g 2 E C F N c Z 2 



+ 



Jkr 



k,q 



g 2 C F /3Z 2 
'3-D 



^ikr 



k 2 + 



nir 



(k 2 + m 2 ) 2 



q 2 + m| 



Ami 



(k 2 + m 2 )[(k + g) 2 + m 2 ]q 2 (k 2 + m 2 ) 2 [(&: + q) 2 + m 2 ];? 2 



+ 



1 



k 2 + 



g 4 C F N c 



(3Z 2 



(k + q) 2 + m| 

„ikr \ 2 

+ c% 5 



(3.18) 



/ k k 2 + m 2 / 

In the terms of 0(<? 4 ), we can immediately set g 2 — > g 2 , Zi,Z 2 —> 1. 

Let us inspect the role that eq. (|3.18p plays when first subtracted in an unresummed form, 
and then added "as is", to eq. (|3.15[) . The first line of eq. (I3.18D accounts for the first two 
lines of eq. (13. 15H . cancelling the infrared divergence like in eq. (I2.13P and fixing 



1 + 



(47T) 



UN* 



L h + 1 



2N { ( 



3 V 



Uf-i 



(3.19) 



The expression in the curly brackets in eq. (|3.18p agrees, for m E — > 0, with the zero-mode 
part of eq. (|3.5p , and through the subtraction-addition step replaces it with a "less" infrared 
sensitive expression; we return to this term presently. Finally, the last term of eq. (|3,18p can 
be written, after subtraction and addition and insertion Z 2 — > 1, as 



-Mr) 



+ 



unresummed 



+ 



9 4 C F N C 



1 exp(— 2m E r) 



8T 2 r 2 



8T 2 r 2 



(3.20) 



rcsummcd 

Now, we observe a problem. After combining eq. (|3,15p with the re-processed version 
of eq. (|3.18p . as outlined above, the large-distance behaviour of ip c is dominated by an 
uncancelled power-law term, 



7rTr>l g 4 C F N c 



1 



+ 



1 



8TV 



9 4 C F N C 



(47T) 



1 



12T 2 r 2 



(3.21) 



(4vr) 2 [ 24T 2 r 2 

Indeed all other terms are exponentially suppressed after the subtraction-addition step, either 
as exp(— 2irTr) or (from eq. (|3.18|) ) as exp(— m E r) or exp(— 2m E r). A term of the type in 
eq. (|3.2ip must be a gauge artifact: physically there is a finite screening length in a non- 
Abelian plasma§| Indeed we will find in sec.[5](cf. eq. (|5.9p ) that in gauge-invariant observables 
the power-law terms do get duly cancelled. 

8 Within perturbation theory a power-law term could in principle also indicate a sensitivity to colour- 
magnetic modes; however, as eq. (|3.20p suggests, the term here has at least partly a colour-electric origin. 
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Despite the issue of eq. (|3.2ip . which we consider to be a serious one from the concep- 
tual point of view, we wish to complete in the remainder of this section our discussion of 
eq. (|3.18p . This leads to another issue, yet this time more physical, being a manifestation of 
the logarithmic sensitivity of the result to the non-perturbative colour-magnetic scale. 

First of all, if we just carry out the integrals in the curly brackets of eq. (|3.18j) literally, the 
result appears to be both infrared and ultraviolet finite, and reads 

-2/^»_ -y2„— mpx i / s3 (~< — „—m E r \ 2 



const. + 



9<C F Z(e- 
AnTr 



+ 



1 fg 2 c F e~ 



4irTr 



+ 



g A C F N c e~ m ^ 



2 - ln(2m E r) - 7e + e 2mEr J Bi(2m B r 



g 4 C F N c exp(-2m E r) 



(4vr) 2 8T 2 r 2 
r»7r/9 2 T g 4 C F N c e- mEr 



+ 0(g 5 



ln(m E r) + 2-ln2-7 E + C 



1 



(3.22) 



+ 0(g 5 ) . (3.23) 



{An) 2 [_ vm E r 
On the last row we displayed the large-distance behaviour. However, the logarithmic depen- 
dence on r (see also ref. |23| ) implies that this term cannot be interpreted as a mass correction 
to the leading order result, but that the correction of 0(g 4 ) overtakes the correction of 0(g 2 ) 
for r ^> ir/g 2 T, rendering the perturbative series out of control. A further issue is that the 
gauge dependent term, the last one within the curly brackets in eq. (|3.18p . does give a fi- 
nite non-zero contribution to eq. (|3.22p . These well-known issues have lead to attempts at 
interpreting the self-energy insertion in a different way [23], and we now discuss these. 

The idea is to treat the 1-loop self-energy as if it were an analytic function of k 2 , as would 
be the case if there were no infrared problems. In this case, making use of the symmetry in 
k — > — k, we can write the self-energy contribution (curly brackets in eq. (I3.18P ) as 



Mr) 



+ g 2 c F (3 
g 2 c F (3 



d 3 k 

(2vr) 3 



,ik-r 



1 



k 2 + ml 



+ 



+ 



Air 2 ir 

g 2 c F 

AirTr 



dk k e 



kr 



(k 2 + m 2 ) 2 
IIo (im E ) + (k- im E )UQ Q ' \im E ) + ... 



1 



2im E 



k 2 + 



exp 



mi 



(k 2 + 



in 



2\2 



m E + 



00 



2m P 



(3.24) 



where we closed the contour in the upper half plane and then resummed the mass correction 
into the leading order term. Inserting the values (cf. eq. (|2.17|) ) 



n oo(^E) 
2m E 

no \im E ) 



g 2 N c T(_i 



87r \e m . 



In 



4m 2 



+ 1 



n E ' 



2inriT, 



0, 



(3.25) 
(3.26) 
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we note that handled this way, the dependence on the gauge parameter £ disappears; however, 
an uncancelled infrared divergence, associated with contributions from the colour-magnetic 
scale, is left over. In fact the term in eq. f|3.25j) equals the ultraviolet matching coefficient 
from the scale m E that plays a role in the approach of ref. |4Q| (with the only difference 
that a subtraction of an unresummed infrared contribution, g 2 N c T/8-K x (l/e uv — 1/ £ ir) = 0) 
transforms l/e m into l/e uv in that case). 

The third possibility is to use the pole mass method but to "regulate" the infrared behaviour 
of the spatial gluon propagators by introducing a "magnetic mass" as a regulator (see, e.g., 
refs. [SUES]). In this way a finite and gauge-independent result is obtained; however, it is 
ambiguous, because the value of the magnetic mass has no well-defined meaning. 

On the non-perturbative level, we expect that the infrared problem is cured by physics at 
the colour-magnetic scale g 2 T/ir; this physics being that of three-dimensional confinement, 
it however cannot be reduced to a simple magnetic mass. Nevertheless, we could still ex- 
pect exponential decay as in eq. (|3.24p . with a non-perturbative version of the correction in 
eq. (|3.25p . that is with the screening mass 



If the mass m E were dictated by the mass of the lightest gauge-invariant state, in the spirit of 
ref. [10] . which in perturbation theory reduces to the structure of eq. (|3.25p . then c E sa 6.9 for 
N c = 3 [H]. This is just a guess, however; the result can as well be gauge-dependent, given 
that the large distance r-dependence of the Coulomb gauge correlator appears in any case to 
be determined by an unphysical power-law contribution, as discussed above. Therefore we 
treat m E as a free parameter for now. 

In fig. we compare the magnitude of the 0(<7 E ) correction with that of the O(g^) term, 
with the former one computed either as in eq. (|3.22p . or as in eqs. (|3.24p - (|3.27p . It appears 
that if we focus on small distances, r < 2/T, then the "unresummed" method, i.e. eq. (|3.22p . 
shows better convergence. Therefore in the following we mostly concentrate on this case 
although, as the discussions above and in fig. 0] show, the unresummed result cannot be 
extrapolated to large distances, r > 2/T. 




(3.27) 
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Figure 4: The 0(g E ) correction to i/'c (the expression with curly brackets in eq. (|3.22p ) over the C(.g E ) 
term (the first r-dependent term in eq. p.22|l ). labelled as the "unresummed ip®"; together with the 
same ratio emerging from the procedure of eqs. (|3.24j) - (|3.27|) . labelled as the "resummed" method. 
At small distances, the unresummed method shows better apparent convergence, particularly if m E is 
small as would be preferred by fig. [5] (where m E > 2T would lead to stronger screening and a larger 
deviation from lattice data); however, the unresummed method must not be used at r 3> 1/T. 



3.5. Summary and comparison with literature 



For moderate distances, r <C 7r/g 2 T, we can write down a "full result" by adding up eqs. (13.150 
and (13321) : 
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In 



1 + e" 



-27rTra; "1 



1 - e 



-2nTrx 



+ 0(g 5 ), 



(3.28) 



where L^,L{ are as defined in eq. (|3.16p and £?i,Li2 in eqs. (|A.20|) . (jA,21jh This result is 
renormalization group invariant up to 0(g 5 ). For large distances, the square bracket term 
on the second line is to be omitted, and the exponential function on the first line is to be 
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Figure 5: The singlet potential in Coulomb gauge from eq. f|3.28[) \f$V\ = - ln(^ c /|7/i P | 2 )], for Nf = 0, 
at T = 3.75Aj^g (T rs 3T c ), at small (left) and large (right) distances. The band corresponds to 
variations of the gauge coupling and m E as explained in the caption of hg. [2] As explained in fig. |4l 
"unresummed" results can only be applied at short distances, "resummed" ones only at large distances. 
The lattice data, labelled by N T , is from ref. [S] (the spatial lattice size was kept fixed at 32 3 ). 



replaced with exp(— m E r), with some non-perturbative m E . The uncancelled power-law term 
on the third line implies that the singlet free energy dies away at large distances slower than 
gauge-invariant correlations. 

In fig. [5] we compare eq. f|3.28|) with Nf = lattice data from ref. [8]. The parameters 
have been fixed as in eq. ()2.22p . and also more elaborately as explained in the caption. We 
observe good agreement between our result and the non-perturbative data, if the unresummed 
form of eq. (|3.22p is used at short distances, and the resummed form of eqs. (|3.24|) - (|3.27|) 
at large distances. (Unfortunately the latter expression involves an unknown parameter, 
m E , so the test is less stringent at large distances.) We have repeated the comparison at 
T ~ 12T C , and the agreement remains good, despite the band becoming narrower (cf. fig. [6]). 
Such a nice agreement for ip c even at T « 3T C is perhaps somewhat surprising, given that 
according to fig. [2] higher-order perturbative corrections to ip P could still to be significant in 
this temperature range. (Formally, ip P can be obtained from the T- and r-dependent part of 
tp c by setting r — > 0, cf. eqs. (|2.9|) . (|3.2p . and it can indeed be observed from fig. [SJleft) that 
some tendency towards a discrepancy starts to form in this limit.) 

It is interesting to compare the present results with those in ref. |42| . where the short- 
distance spatial correlators related to gauge-invariant scalar and pseudoscalar densities were 
measured. The authors observed stronger correlations than indicated by the leading-order 



18 



N, = 0, T= 12 T 

f ' C 



N, = 0, T = 12 T 

I ' C 



> 

03. 




> 

on. 



-2.0 



0.00 



-0.02 



-0.04 : 



-0.06 




2-loop g , resummed, m E = 1.5T 

2 ^ 

1-loop g , resummed, m E = 1.5T 



1.2 1.4 1.6 1.8 2.0 2.2 2.4 
rT 



Figure 6: Like fig. [5] but for T = 15A^g (T w 12T C ). The parameter m E has been scaled down 
roughly by an amount corresponding to the running of the gauge coupling [21J . 



perturbative predictions (in the language of the static potential, this would correspond to less 
screening). It would be interesting to see whether pushing the perturbative determination of 
those observables to the same level as here would allow to make the match as good as that 
in fig. [5j 



4. Singlet free energy in a general covariant gauge 



The singlet free energy in a general covariant gauge is defined in analogy with eq. (|3.ip . 



^^(Tr^P^-. 



(4.1) 



The graphs are the same as in fig. 02 and so is the general unresummed result of eq. (13.2 
The difference is that for the gauge field propagator we now insert 



(4.2) 
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while the self-energy has the form 

n 00 (o,k) 
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where -?T = (0,k). Expanding in k 2 following the philosophy of eq. (I3.14p . the self-energy 
behaves as 
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where terms of 0(k 4 /T 2 ) are ultraviolet finite. The fermionic part here has a structure 
familiar from eq. (13.151) . 

While in principle we could go on as before, working out the contribution of 0(k i /T 2 ) in 
detail, it seems to us that the result is not particularly interesting. This can be seen already 
at short distances, mT <C 1, focussing on the divergences, as we will now do. 

Listing only the 1/e-poles of dimensional regularization, the renormalization of the bare 
gauge coupling (cf. eq. (|2.3h ) from the first graph of fig. [3] yields the contribution 
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while the fourth and fifth graphs are ultraviolet finite. The sixth graph yields 
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and the seventh graph, inserting the expansion in eq. (|4.4p . amounts to 
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Adding up, we observe that divergences do not cancel, unlike in Coulomb gauge, but that 
the results sum up to 
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In principle we could "renormalize" the result by writing 
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however the form is not that of known endpoint divergences (see, e.g., ref. |43| ) which are 
additive, rather than "correcting" the renormalized g 2 ; in addition, the finite parts of the 
potential still remain gauge dependent, even at short distances mT <C 1, because of the 0(1) 
term in eq. (|4,10p . Therefore the result does, in general, not match the zero-temperature 
potential. The situation is analogous to that in the next section, where it will be discussed 
somewhat more explicitly (cf. eqs. (15 .7p . (|5,8p ). 

Let us point out that, compared with the conventional zero-temperature static potential 
where divergences do cancel [38], the only difference is with the second and third topologies 
in fig. [3l At zero temperature, the time direction is infinite and the Wilson loop has a finite 
extent within it; at finite temperature, the time direction is finite and the Polyakov lines 
wrap all the way around the time direction. This turns out to lead to a difference for these 
particular topologies. 

5. Cyclic Wilson loop 
5.1. Basic setup 

With the objects in eqs. (|2.ip and (|2.2p . the cyclic Wilson loop is defined as 



where we made use of periodic boundary conditions (i.e. Wp = Wo). In an Abelian theory 
V>w would agree with ip c of eq. (|3.ip . whereas in the non- Abelian case ip w can be viewed as 
a gauge-invariant "completion" of i/j c - 

In general, the insertion of the spacelike Wilson lines into the observable leads to a huge 
proliferation of graphs compared with those in fig. [3j On the other hand, the fact that both 
W and W f appear, i.e. that the (untraced) Polyakov loops are really connected by a Wilson 
line in the adjoint representation, also means that there are many cancellations. Even so, 
quite a number of diagrams remains. To simplify the task somewhat, we have restricted to 
special gauges in this case, namely to those where the gluon propagator has no components 
mixing time and space indices, i.e. Goi{Q) = 0; this class includes, in particular, the Coulomb 
and the Feynman gauges. With these provisions, only the graphs in fig. [7] give an additional 



^w(r) = ^(TV[P r W P o t W t ]) 
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Figure 7: The additional graphs relevant for the cyclic Wilson loop, supplementing those in fig. G2 
in the class of gauges where the gluon propagator has no components mixing time and space indices, 
i.e. Goi(P) — 0. Time runs vertically and space horizontally. 



contribution. Choosing for convenience the coordinates so that r points in the z-direction, 
the general (unresummed) result can be written as 



/e- G „„(0, k )/(e--l)E 

1 Qn 



G zz (q n ,q) 



+ 



k,q 



e -iq r _ g ik r 
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k z + q z 



(5.2) 



where ip c refers to the general result in eq. (|3.2|) . It can be verified, by taking the Coulomb 
and Feynman gauge results for ip c from sees. El HI respectively, as well as the corresponding 
propagators from eqs. (|3.3p . (|4,2p . that the expression in eq. (|5.2p is indeed gauge independent. 
(In fact, we will demonstrate the gauge independence of the EQCD part of the expression 
around eq. (|5.12p below.) It is also trivial to see that formally the new contribution in eq. (|5,2p 
vanishes for r = 0, like the structures of eq. (|3.2p . 

Inserting the Coulomb gauge propagator from eq. (|3.3p and carrying out some changes of 
integration variables, we can write the unresummed result as 
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(5.3) 



Now it is understood that the Coulomb gauge ip c from eq. (|3.5p is to be inserted. 



5.2. Short-distance limit 

Let us start by inspecting the new terms in eq. (|5.3p at short distances, in analogy with 
sec. 13.21 The latter row in eq. (j5.3j) only contains the Matsubara zero mode and, compared 
with the contribution from the sum, leads for dimensional reasons to a contribution sup- 
pressed by rT at small distances (once (3 is factored out). The zero-temperature integrals 
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corresponding to the structures on the first row can be carried out, 
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The divergences add up; given that t/> c was finite after charge renormalization, there is noth- 
ing more available that could cancel them. We may, nevertheless, represent their effect by 
introducing a "renormalization factor" Z w in analogy with eqs. (|4.10p . (|4.1ip . 

ikr 



In 



Z w g z C F p 



k 2 



+ 0(g 4 



(5.6) 
(5.7) 
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but this does not have the form normally related to the cusps that appear in the observable [26] 
(the corresponding divergences are additive, rather than "correcting" the renormalized g 2 ). 
In addition, there is an ambiguity concerning whether the Fourier transform in eq. (|5.6p is 
to be understood in 3 — 2e or in 3 dimensions. Preferring therefore to write down a bare 
expression, we can summarize the short-distance behaviour as 
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where we made use of footnote [7] in order to rephrase all the logarithmic dependence in lc- 
space. The expression in eq. (|5.8p is gauge-independent; however, it does not go over into the 
zero-temperature potential at short distances. 



5.3. Hard-mode contribution 

In order to determine the behaviour of V\v a t larger distances, r ~ 1/vrT, we treat the new 
terms in eq. (|5,3p according to the strategy around eq. ()3.14p . Two of the sum- integrals are 
by now familiar, and given by eqs. (|A.7p . (|A.15|) : the new one is given by eq. (|A.19p . and we 

get 
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81n(27rTr) 



(5.10) 



Here we already replaced the sum over the exponential integral through a simple integral 
representation, cf. eq. (|A.22|) . 

At short distances, nTr <C 1, the complicated square bracket expression in eq. (|5.9p goes 
over into 

g A C F N c 

(4vr) 3 Tr 

Combining with the logarithmic term from the first row (with inserted from eq. (|3,16p ). 
we reproduce the result of eq. (|5.8p . At large distances, on the other hand, the power-law 
term in eq. (|5.9p cancels against that in eq. (|3.2ip . Therefore, apart from a disconnected part 
like in eq. (|2.5p . ip w (r) is exponentially suppressed at large distances. 

5.4. Soft-mode contribution 

We denote the EQCD contribution to ip w by if^. Within EQCD, the new graphs in fig. [7] 
amount to an evaluation of the expectation value (see also ref. |35j ) 

Kir) = r c (r) + £^(y*Tr [A (r)W A (0)W^] - Tr [A (r)A (0)]) + 0(g 5 ) , (5.11) 

where the subtraction corresponds to the part already included in We have introduced 
another "iJ-factor", this time denoted by yo, related to the fact that the spatial Wilson lines 
within EQCD might differ in normalization from those in QCD. The new contributions we 
are interested in are (D(g 4 ), so that Z\ can be set to unity and g 2 can be set to g 2 in the 
evaluation of eq. (|5,lip , A straightforward computation leads to 

^(r) = r c {r) + g 2 c F i5(yl 
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(5.12) 

where we have for completeness kept a general gauge parameter. It is easy to check now 
that the ^-dependent part cancels against the self-energy contribution to ip® from the 1st 
row of eq. (|3.24j) . with the self-energy inserted from eq. (|2.17p . so the result is indeed gauge 
independent. 

On the other hand, if we take the Coulomb gauge result as a starting point, then all terms 
in eq. ()5. 12[) need to be kept, with the value £ = —1. The term on the 2nd row of eq. ()5.12p 
is factorized and contains the integral 

r e iQzr _i r ( 1 r 1 1 

/ 2 2 =-— — + 2 hu>)+ 7E -l , (5.13) 
J q qj q l 8vr [ e uv L J J 
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multiplied by the usual g A C F N c J k e ik r /(A; 2 + 



g A C F N c ex.p(-m E r)/4:irr + 0(e) (cf. 



eq. (|5.20p ). It is important to stress that, as underlined by the notation, the divergence in 
eq. (|5.13p has an ultraviolet origin. This can be seen, for instance, by regulating the ^-integral 
by taking a principal value and the gj_-integral by introducing a mass- like regulator: 
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The ^-independent part of the square brackets is multiplied by 
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(5.15) 



whereby we recover the divergent part of eq. (j5. 13|) . This divergence is essentially the self- 
energy correction related to an adjoint Wilson line |40j . 

As far as the other terms in eq. ()5. 12j) go, the last one simply yields 
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The remaining (middle) term requires a bit more work; some intermediate steps are given in 
eqs. (|A,29p - (|A.3ip . The final result can be written as 
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It may be noted that the large-distance behaviour of eq. (|5.18p largely cancels against that 
in eq. (j332|) . 

To summarize, the additional contribution to is composed of the first new term in 
eq. (|5.12p . together with the results from eqs. (|5.13p . ()5.16p . and (|5.18p . At small distances 
the "old" part V>c could be taken from its direct evaluation, eq. (|3.22p . while at large distances 
the resummed form of eqs. (|3.24p - (|3.27p is preferable. 

It remains to fix in the first new term of eq. (|5.12p . In order to match the behaviour 
on the 2nd row of eq. (|5.9p in the subtraction-addition step, we need to choose 
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It can be seen, however, that there is a potential ambiguity from terms of the type 0(l/e) x 
0(e) that this multiplies. Here the 0(e)-terms come from the massive or massless leading- 
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order potential, 
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where we inserted 1 = /i~ 2e [l + e(ln ^- + 7 E )] in order to fix the dimensions. It is not clear to 
us whether the 0(l/e) x 0(e) terms from here can have physical significance. 

In any case, after fixing y$, the terms can be added up. The complete result is not 
particularly transparent, and may be ambiguous as just discussed, so we do not write it down 
explicitly; it suffices to say that the sum can be expressed as 
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The function Q Dn , in which the complications are hidden, does have a simple expression in 
certain limits, however, and these will be discussed in the next section. 



5.5. Summary and comparison with literature 

Adding up the relevant parts of eqs. (|3.28j) . (|5.9p and (|5.22p . we finally get 
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At small distances, m E r <C 1, the dominant term of the coefficient function <5 DR , defined in 
eq. (I5T221) . is 
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Though we have not carried out the computation, we could expect that in lattice regulariza- 
tion the corresponding structure goes over into 
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where g 2 is a suitably defined renormalized coupling. Therefore, it would appear that Q UI 
diverges in the continuum limit. On the other hand, at large distances, m E r ^> 1, the 
dominant term of the coefficient function is 
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where possible logarithms of m E r are also included in 0(1) (whether such logarithms appear 
is related to the ambiguities mentioned after eqs. Apart from the said 

logarithms, eq. (I5.26D can be accounted for by a mass correction, 

exp(— m E r) — > exp(— m E r) , (5.27) 

where m E contains a logarithmic ultraviolet divergence: 

g 2 N c T( 1 p? 
h in — 

-UV m E 

A naive transliteration to lattice yields 



m E , DR = m E + ^-^ + 1*^ + 0(1)). (5.28) 



m E , lat = m E + fin — + (9(1)) . (5.29) 

87r \ a z m E / 

So, Debye screening grows logarithmically as the continuum limit is approached, much like 
in the numerical study of ref. [H], and the exponential function exp(— m E lat r) decreases. 

To summarize, we believe that the leading term, ~ £/ lat exp(— m E r), is a very ultraviolet 
sensitive and pathological function of the lattice spacing; at large distances, it appears to ex- 
trapolate towards zero in the continuum limit, while at short distances it appears to explode. 
In fact it might look somewhat like a delta-function. 

Concerning whether the same happens also in the other terms, in particular in the one 
with exp(— 2m E r), a higher order computation would be required to see any perturbative 
indications. We find it conceivable, though, that some of the terms could also remain finite, 
representing a coupling to the gauge- invariant channel of the traced Polyakov loop, 

Mr) = ^(Tr[P r ]Tr[P^}). (5.30) 

This behaves as 123 



f ip T (r) \ g 4 C F exp(-2m E r) 5 
HWF J {^?N C 4T 2 r 2 + ° {9 } (5 - 31) 

at the order of our computation; the argument of the large-distance exponential fall-off has 
recently been determined at next-to- leading order in ref. [30 1 . (A more precise analysis of 
this correlator at short distances has been undertaken in ref. [28J.) 



9 No logarithms of r appear if the soft contributions to ipc are treated as in eq. (|3.24[) and no terms of the 
type O(lfe) x 0(e) are included, i.e., if we just sum together eq. ()5.13p and the large-r limit of eq. (|5.18p . 
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Figure 8: Two sketches for the potential extracted from the cyclic Wilson loop, /3V cyc \[ c = 
— ln^w/l^pl 2 ), after subtracting the strongly cutoff-dependent part A oc Q r exp(— m E r) from 
eq. ([5723]) . The plot is for N { = and T = 3.75Ams i T ~ 3T C ). The band corresponds to varia- 
tions of the gauge coupling and m E as explained in the caption of fig. [U Comparing with fig. [SJleft), 
we see significantly less correlation (stronger screening). 

A numerical sketch of eq. (|5.23p . obtained by omitting the regularization dependent leading 
term altogether, and including two versions illustrating the importance of the remaining 
EQCD contribution, is given in fig. [8) The main "prediction" we can make is that in the 
continuum limit the result should be much closer to zero (i.e. much more strongly screened) 
than for the observable in fig. [3j Unfortunately we cannot compare this prediction with the 
lattice data from ref. [13], because the spatial Wilson lines were smeared in that study. On the 
other hand, a direct lattice measurement of the cyclic Wilson loop does appear to reproduce 
both the strong lattice spacing dependence as well as the small continuum extrapolation that 
are suggested by our analysis [44|. 

6. Conclusions 

The purpose of this paper has been to compute, within the weak-coupling expansion in con- 
tinuum, a number of Polyakov loop related observables that have been popular in the context 
of phenomenological studies related to the fate of heavy quarkonium at high temperatures. 
Our main findings can be summarized as follows. 

At large distances, many of the observables considered go over to a gauge- independent 
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disconnected contribution, the square of the expectation value of a single Polyakov loop. We 
have recomputed this quantity at 0{a%) and found a result different from the classic one in 
the literature (we believe that ours is the correct one). Numerically the difference is very 
small and does not improve on the match to lattice data; for that, inclusion of higher order 
terms in the weak-coupling series would probably be needed. 

The singlet free energy in Coulomb gauge, defined by eq. (|3.1|) and denoted by ip c , is ultra- 
violet finite without any additional renormalization factors in dimensional regularization at 
order O(a^), and agrees at small distances with the gauge-independent zero-temperature po- 
tential [3S]. At large distances, however, its behaviour is not determined by Debye screening, 
but by a power-law contribution, which in our view is related to gauge fixingl^j In principle 
this means that the Coulomb gauge potential could be more binding than would physically be 
expected. On the other hand, at intermediate distances the power-law term has no dramatic 
effect. In any case, our result does agree with lattice data surprisingly well (cf. fig. [5]), even at 
low temperatures; a higher-order perturbative computation could perhaps tell whether this is 
an accident or a robust characteristic of thermal correlations at non-zero spatial separations. 

The singlet free energy in a general covariant gauge contains ultraviolet divergences and 
remains gauge dependent even at short distances, where it in general does not match the 
zero-temperature potential. 

The cyclic Wilson loop is, by construction, extracted from a gauge-invariant quantity, and 
indeed the gauge-fixing related power-law term affecting ip c cancels out. However, the be- 
haviour of this observable does not match the zero-temperature potential at short distances. 
Moreover, at large distances it is not sensitive to physically Debye-screened one-gluon ex- 
change, but rather develops a logarithmically ultraviolet divergent screening mass. Neverthe- 
less, there could remain some non-trivial higher-order r-dependence in the continuum limit, 
either with a screening mass ~ 2m E or with ~ 27rT. 

Our findings underline the pitfalls that exist in generic potential model studies of quarko- 
nium physics at finite temperatures. On the other hand, it can be argued that parametrically 
only the range rT < 1 is relevant for quarkonium melting [46 j , so problems that are related 
to the asymptotic large-distance behaviour may not be that important from the pragmatic 
point of view. With this philosophy the Coulomb gauge singlet free energy seems like the 
least misleading quantity to use, if it is assigned the meaning of the real part of the proper 
potential [6]. At the same time, it appears that in this range perturbation theory alone 
may yield a fairly good description of the system, so that one could just as well resort to 
gauge independent perturbative quarkonium spectral function determinations such as those 
in ref. [37] . 



It is interesting to note, though, that ours is not the first study suggesting the existence of power-law 
terms in heavy-quark related observables, cf. e.g. ref. [45] . 
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Appendix A. Sum-integrals and integrals 

In this appendix we list, for completeness, the results for the four-dimensional sum-integrals 
and the three-dimensional integrals that play a role in our computation. The theory is 
regularized dimensionally, with D = 4 — 2e, and terms of 0(e) are omitted. The dependence 
on the scale parameter, fl, appears in most cases through logarithmic factors which are 
denoted by L\>,Lf and have been defined in eq. (|3. 16j) . A prime in the sum-integral symbol 
indicates that the Matsubara zero mode is left out; when the summation- integration variable 
is in curly brackets, the Matsubara frequencies are fermionic. The four- vector denoted by K 
is by definition purely spacelike, 

K = (0,k). (A.l) 
With this notation, the following sum-integrals can be determined (Q 2 = + q 2 ): 
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is an exponential integral (£1(2:) = — Ei(— z) for z > 0), and 
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is a dilogarithm. Where not shown explicitly, the fermionic sum-integrals are obtained from 
the corresponding bosonic ones by replacing Lb — > Lf. 

In eqs. (|A. 13^ — (|A.19|) . various sums over the exponential integral appear. All of these can 
be transformed into an integral representation: 
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Within EQCD, the following integrals can be worked out: 
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Finally, we consider the EQCD-integral defined in eq. (|5.17p . In fact, it is useful to re-start 
with manifestly infrared safe integration variables like in eq. (|5.2p . and then the integral can 
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be expressed as 
X = 
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Here we first substituted q — > — k — q; the numerator vanishes for q z = so the integrand 
is regular. Subsequently we changed integration variables as k — > k — q and q — > — q in the 
second term, introducing a principal value to regulate the potential divergence at q z = 0. 
The second term on the last line can be written as 
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where in the last step the integration could be carried out with the help of a Feynman 
parameter. Combining with the first term of the last line of eq. ()A.29[) . we get 
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167r 2 m E r 

where we proceeded as in connection with eq. (|3.22p . 
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